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Abstract: Let B H = {B H (t),t E M. N } be an (TV, d) -fractional Brownian sheet with 
Hurst index H = (Hi,... ,Hpj) E (0,1)^. The main objective of the present paper is 
to study the Hausdorff dimension of the image sets B H (F + t), Fa R N and t E R N , 
in the dimension case d < + • • • + . Following the seminal work of Kaufman 

[9], we establish uniform dimensional properties on B H , answering questions raised by 
Khoshnevisan et al. [11] and Wu and Xiao [17]. 

For the purpose of this work, we introduce a refinement of the sectorial local- 
nondeterminism property which can be of independent interest to the study of other 
fine properties of fractional Brownian sheets. 
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1. Introduction 

In the last thirty years, several extensions of the well-known fractional Brownian motion 
(fBm) introduced by Mandelbrot and Van Ness [12] have emerged in the Gaussian random 
fields literature. Two major classes of multiparameter processes have been defined: Levy’s 
IV-parameter (fractional) Brownian motion and fractional Brownian sheets. The first one 
is an isotropic process known to be be locally non-deterministic (LND, see Pitt [15] for a 
more complete reference), self-similar and with stationary increments. As a consequence, 
the geometry and fine properties of the IV-parameter fractional Brownian motion have 
been extensively documented by extending the classic techniques developed in the literature 
related to the fractal geometry of the one-dimensional fractional Brownian motion. 

On the another hand, the understanding of fractional Brownian sheets (fBs) introduced 
by Kamont [8] has been proved to be more challenging and technical as the former does 
not satisfy the classic LND property. Recall that an ( N , d)-fractional Brownian sheet B H = 
{B H (t),t e with Hurst index H = (Hi ,..., Hjy) e (0,1)^ is defined as a centered 
Gaussian process with independent and identically distributed components whose covariance 
is given by 


N 

E[B*( a )B*(t)] = n[N 2 ^ + \U\ 2H ‘ - k - U\ 
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Note that similarly to fBm, it also admits an integral representation with respect to the 
Brownian sheet W: 


B, 


H 




dW„ 


( 1 . 2 ) 


In the case H\ = ■ ■ ■ = Hn = \, we obtain the well-known Brownian sheet. 

Anisotropic Gaussian random fields such as IBs have raised an increasing interest in 
recent years as they appear naturally in the study of stochastic partial differential equations 
(SPDEs) and Markov processes [14]. In a more applied perspective, several phenomena in 
image processing, hydrology and spatial statistics [5, 2] are intrinsically anisotropic, and 
thus, required the introduction of such theoretical models. 

The study of distributional properties of fractional Brownian sheets have been consid¬ 
erably eased with the introduction by Khoshnevisan and Xiao [10] of the so-called secto¬ 
rial local-nondeterminism property. Namely, the latter states that for any u, v, t 1 ,..., t n e 
[£,+ 00 )^, 

N 

'Var(B H (u)\B H (t 1 ),...,B H (t n ))^c 0 'y y min | U# — tj \ 2Hi (1.3) 

1=1 

where the constant cq > 0 only depends on e. Firstly introduced on the Brownian sheet, it 
has been then extended to general fractional Brownian sheets by Wu and Xiao [16]. 

This sectoral LND has been the cornerstone to the study of multiple geometrical proper¬ 
ties of (fractional) Brownian sheets, allowing to adapt classic techniques used on multiparam¬ 
eter fractional Brownian motion to this class of processes. More precisely, the distributional 
properties of the local time and level sets have been investigated by Khoshnevisan and 
Xiao [10], Ayache et al. [1], extending earlier works by Dalang and Walsh [4], Xiao and 
Zhang [21]. The fractal geometry of image sets B H (F) has also been extensively studied 
in [11, 10, 18, 16]. Note that as pointed out by Xiao [20], due to the anisotropic nature of 
fractional Brownian sheets, it is usually convenient to study geometrical properties using 
the following anisotropic metric p : 

N 

Vs, t e R n ; p(s,t) = ^e-te\ H ‘. (1-4) 

i=i 

As previously outlined, we aim in this work to investigate uniform dimensional properties 
of image sets B H (F) of fractional Brownian sheets. The high dimension case d ^ + 

-b -gt- has been thoroughly discussed by Khoshnevisan et al. [11], Wu and Xiao [18] who 

obtained the following result: with probability one, 

for every Borel set F <= R^; dim H I?' ff (F) = dim^F, (1.5) 

where dim([ designates the Hausdorff dimension with respect to the anisotropic metric p. 
Note that a similar result exists on the multiparameter fractional Brownian motion when 
ff ( see the work of Monrad and Pitt [13]). 

When d < ^ + ''' + the previous uniform result does not hold any more. For 
instance, it is obviously false if one considers the level set F = W -1 (0), where W is a 


Paul Balanga/Image sets of fractional Brownian sheets 


3 


Brownian sheet: Khoshnevisan and Xiao [10] have proved that F has positive Hausdorff 
dimension when d < 2 N, whereas we clearly have dim H W(F) = 0. Nevertheless, following 
the ideas developed by Kaufman [9], one may hope to establish a weaker uniform Hausdorff 
dimension result. Indeed, the former has proved that a one-dimensional Brownian motion 
satisfies a slightly weaker property: with probability one, for every Borel set FcR, 

dim H -B(F + t) = 2dim H -F for almost all let. (1-6) 

Khoshnevisan et al. [11], Wu and Xiao [18] have investigated the extension of this property 
to the (TV, 1) Brownian sheet and fractional Brownian sheets satisfying H^d ^ 1 (assuming 
that Hi < • • • < Hjh). Even though it may seem to natural that the former result would 

hold for any fBs such that d < H-1- , this question was left opened in the previous 

works as the authors observed that techniques based on sectorial LND do not seem to scale 
well the general case (on contrary to the (TV, d)-fractional Brownian motion considered by 
Wu and Xiao [16]). 

Consequently, the main purpose of this work is to close the gap between the statement 
of Wu and Xiao [18] and the uniform case presented in Equation (1.5). More precisely, we 
prove in Section 3 the following two uniform results on the geometry of fractional Brownian 
sheets. 

Theorem 1. Let B H be a fractional Brownian sheet and suppose d < X!^=i TT • Then, with 
probability 1, for every Borel set F ( 0 , 00 )^, 

dim h B h (F + t) = minjd,dim^E 1 } for almost all t e . (1-7) 

In addition, we also extend the result of Wu and Xiao [18] related to Lebesgue measure 
of image sets. 

Theorem 2. Let B H be a fractional Brownian sheet and suppose d < Xirli ~rr t - Then, with 
probability 1, for every Borel set F <z (0,00)^ such that dim^E > d, 

A d(B H (F + t)) > 0 for almost all t e , 

where A d denotes the Lebesgue measure on M. d . 

The proof of the two previous results rely on the introduction in Proposition 3 (Section 2) 
of an anisotropic local non-determinism property different from sectorial LND. The former 
then allows to adapt the seminal methods of Kaufman [9] to fractional Brownian sheets 
satisfying d < TT t - N°t e that we hope that this anisotropic LND property can also be 
of independent interest to the study of remaining open questions on the fractal geometry of 
fractional Brownian sheets and more general anisotropic Gaussian random fields. 

2. Anisotropic local nondeterminism 

The local nondeterminism property has historically been introduced by Berman [3] in the 
study of local times of Gaussian processes. Since then, it has been widely and successfully 
used to obtain multiple fine sample paths properties of Gaussian processes, including small 
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balls probabilities, level sets and Hausdorff dimension of graphs and image sets. We refer to 
the surveys of Geman and Horowitz [6], Xiao [19] for a more precise overview on the subject. 

As previously outlined, the Brownian sheet, and thereby fractional Brownian sheets, are 
known to be non locally non-deterministic and Khoshnevisan and Xiao [10] have introduced 
the sectorial local nondeterminism property in order to still being able to investigate distri¬ 
butional properties of this class of processes. The simplest form has been presented in the 
introduction, Equation (1.3). In order to investigate uniform dimension of image sets, one 
needs an analogue of the former on increments. Namely, Wu and Xiao [18] have proved that 
any fractional Brownian sheet satisfies for every s, f, s 1 ,..., s n e [e, +oo), 

Var (B H (s)-B H (t) \ B H (s 1 ),..., B H (s n )) 

N r 

^ ci V min-! min \se — s J f \ 2H ‘+ min \te — si\ 2He ,\se — te\ 2Hi >, (2.1) 

(l^jsSn 1 lsy^n 1 J 

where the constant Ci > 0 only depends on e. 

Nevertheless, it appears in the work of Khoshnevisan et al. [11] that the previous sectorial 
local nondeterminism property is not sufficiently fine to extend the result (1.6) to any 
fractional Brownian sheet satisfying d < JT, • Consequently, we present in the following 

proposition a refinement of the former. 

Proposition 3. Suppose B H is a fractional Brownian sheet and e > 0. Then, there exists 
a constant c o > 0 such that for every t, s, s 1 ,..., s n e [e, 1)^. 

Var(B H (t)-B ff (s) \ B H (s 1 ),..., B H (s n )) > c 0 p(s,t) 2 ■ , (2.2) 

e =x 

where for any i e {1,..., N}, we define 

rn := min \st — sl\ + min \tg — sl\. 

1 1 1 * 

Proof. We aim to prove a property slightly stronger than Equation (2.2). Namely, for every 
t, s, s 1 ,..., s n e [e, l) w , 

Var (B H (t)-B H {s) \ B H (s 1 ),..., B H (s n )) 

N N r , 

^ Co 2 min {rj H \\st - ti\ 2He } + c 0 r\ Hl • | - ti\ 2Hi i. (2.3) 

1=1 t= 1 ' 

We easily observe that since min {r 2Hk , |sfc —tk\ 2Hk } ^ r 2Hk • |sfc — tk\ 2Hk , the former clearly 
induces Inequality (2.2). In addition, in the two components appearing in Equation (2.3), the 
first one is a clear consequence of the sectorial LND property (2.1), since = minj^ n |s^ — 
s\ | + min^ n |^ - sj |. 

Hence, we may focus on the second part, and, set k e {1,..., N} and i A k. Without 
any loss of generality, we may assume that t, — Sj ^ 0 (unless, simply permute t. and s). If 
|sfc — tk\ ^ Tk, we simply observe that 

mm{r 2 k Hk ,\s k - t k \ 2Hk } = r 2 k Hk ^ r 2 k Hk ■ |s* - ti\ 2Hi . 
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The combination of the previous remark and the sectorial LND property yield the expected 
inequality. Therefore, we may assume in the sequel that |sfc —tk\ < r k - 

Owing to the integral representation (1.2) of fractional Brownian sheets, we know that 

Var (B H (t)-B H (s) \ B H (s 1 ),..., B H (s n )) 


= inf E 

aeR™ 


B H (t)~B H (s)-Y 1 a j B H (s^y 

j =i 

inf f f K(u,t, H) — K(u, s, H) — V] ajK(u, s 3 , H)\ du. 
aeR" Jrjv V “I ) 

}. The previous expression of the 


where K(u,t,H) := YleLiifa ~ u ^)+ t ^ ~ (~ u e) 
conditional variance can be lower bounded by 


o =i 
H e - l/2-i 


r , N N n N 

f fl(^- u e)+ e ~ 1/2 - “ u f)+ e ~ 1/2 - Yj a 0 fl(4 - u z) 

ae Jr" V ^ =1 i=1 j =1 t=1 


H e - 1/2 


dit. 


In order to obtain a uniform lower bound of the previous expression, the main idea is to 
exhibit an element h e L 2 which is orthogonal to the family of functions D[^ =1 (s^ — iq) + f 1 ^ 2 . 
Hence, let us define 

h{u) = h k (u k ) ■ l[ s . it .](Uj) • Y[ ![(),£] («<)• 

C^i,k 

where h k (u k ) = ( u k - t k + r k ) 1 / 2 ~ Hk + (u k -t k - r k )+ 2 ~ Hk - 2 (u k - t k )+ 2 ~ Hk . Note that 
the support of the function h is included in the set \t k — r k , go) x x [0,£] lV ~ 2 (up 

to a permutation of variables). In addition, when H k = 1/2, h k simply corresponds to 
the difference 1 [t k - rk ,t k ) ~ 1[ t k ,t k +r k )■ Let us prove h is well-designed for our purpose by 
evaluating the scalar product, for any fixed j e {1 ,,n}: 

N 




Hi- 1/2 \ _ 


e=i 


' f (a 

J R-i- 


Uk) 


H k - 1/2 


h k {u k )du k , 


where c corresponds to the integration over variables ue, £ =£ k. We need to distinguish two 
different cases, depending on the value of s 3 k . 

1. If s 3 k < tk — r k , (s 3 k — u k ) H + k 1,2 and h k have disjoint supports (respectively [0, s^] and 
[t k — rfc,oo)), and therefore, the inner product is clearly equal to zero. 


2. If s 3 k >t k - r k , 


N 


t=\ 


Hi- 1/2 


l 




(4 -u k ) Hk 1/2 • h k (u k ) dufc. 


Let a denotes either 0, — r k or r k ■ Then, 

(s 3 k - u k ) Hk ~ 1/2 ■ (u k - t k - a) 1,2 ~ Hk d u k 


1 


[tfc + a,sj] 


= (4 ~t k — a) ■ f v k k 1/2 (1 - v k ) 1/2 Hk dv k , 

J[o,i] 
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using the change of variable v & = (s J k — Uk)/(s k — — cl). Consequently, 


N 


h, m - u ^+‘ 1/2 ) = c \ “ t k - r k) + (si -t k + r k ) - 2(s J k - t k ) [ = 0. 


l= l 


The function h is orthogonal to any — u t)+ e ^ and therefore to the linear space 

spanned by the previous collection. As a consequence, 

Var (B H (t)-B H (s) \ B H (s 1 ),..., B H (s n )) 

N 


> 


1 


N \ 2 

Hi-1/2 y~\, \ Hi —1/2 


1 


h,Y[(te ~ u f)+ / ~YV- Si ~ 

L 2 \ 1=1 1=1 

N 


h ’ Yl^ te ~ u ^+ e 1/2 ^ ’ 

/=i / 


L 2 \ 1=1 


since the support of h does not intersect [0, s] due to the component lr Si;t .](uj) in the former. 


Let us first estimate the norm ||/i||^ 2 : 


|| 2 = £ N 2 ■ |Sj - ti| 


f \(u k - t k + rk) 1 / 2 Hk + (u fc - t k - rk) 1 / 2 Hk - 2(u k - 4) + 2 Hk 1 d u k 
Jr 1 ; 


= e N - 2 -\s^U\-r 2 - 2Hk 


dv k 


= Co |s» 2Hk . 


The previous integral is finite since 1 — 2 H k > —1 and v/ 2 Hk + (v k — 2)/ 2 Hk — 2(v k — 


1) 


1/2 -H k 


-co V k 


(h, Y\(ti - Ut)/ e 1/2 \ = Y[ f (te-ue) Hi 1/2 d uA (ti^Ui) Hi 1/2 d m 

\ p —i / Jst 


-1-2 H k 


. On the other hand, the inner product is equal to 


i,k 
rtk 


I (tk - u k ) Hk 1/2 (u k - tk + r k ) 1/2 Hk 
Jtk—rk 


d U k 


> c 2 I Si - ti\ Hi+1/2 ■ r k . 


still using a similar change of variables and observing that c 2 > 0 only depends on e. Hence, 
we eventually obtain 

Var (B H (t) -B H (s) \ B H (s 1 ),..., B H (s n )) ^ c 3 r 2 k Hk ■ \ Si - U\ 2H \ 

where the constant c 3 > 0 only depends on N, e and H. This last inequality then clearly 
leads to the second term in Equation (2.3). □ 

Remark 1 . We may note that the local nondeterminism property presented in Proposition 3 
is not sensu stricto an extension of the sectorial LND property (2.1). Indeed, one can simply 
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Figure 1: Example of conditional variance Var (B(t) — B(s) | B(s 1 ), B(s 2 )) 


observe that if the two terms p{s,t ) 2 and r i H ‘ are °f same order, then the sectorial 
LND bound (2.1) is tighter. 

On the other hand, Figure 1 illustrates the typical case where the anisotropic LND 
property (2.2) provides a better estimate than the classic sectorial LND. Namely, if W is a 
two dimensional Brownian sheet, we easily observe that the bound given by Equation (2.1) 
on Var(W(t) — W(s) | W(s 1 ),W(s 2 )) is zero. On the other hand, Proposition 3 gives an 
optimal lower bound, proportional to the variance of the term W(t) — W ( s ) — W ( s 2 ) + W{s 1 ) 
(informally equal to Var(VF(“grey area”))). This improvement corresponding to some specific 
geometrical configurations will be the cornerstone in the proofs of Theorems 1 and 2. 

Remark 2. The calculus presented in the proof of Proposition 3 offers an alternative way 
to prove the sectorial LND property. Wu and Xiao [18], and originally Kahane [7], used 
estimates on the Fourier representation of fractional Brownian sheets to obtain the lower 
bound, whereas our proof is based on the classic time integral representation (1.2). 

3. Weak uniform Hausdorff dimension of image sets 

Based on the refinement obtained in the previous section, we now extend the weak uniform 
Hausdorff results presented by Kaufman [9], Khoshnevisan et al. [11] and Wu and Xiao [18]. 
The structure of the proof of Theorem 1 follows the ideas initially described by Kaufman 
[9] and relies mainly on the estimate obtained in the following lemma. 

Lemma 3.1. Define 


I{x, y, R) = l [ _ 1A] (R-\\B H (x + t) - B H (y + t)\\)dt, 


Then, for all R > 0, x, y £ [s, l] jV and integers p ^ 1, 


E[(l(x,y,R)) P ] ^ ^(p\) N R- d Pp(x,y)- d P, 


(3.1) 


where the constant Cq only depends on e. 
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Proof. Since B^ ,..., B^ are independent copies of an (N, l)-fractional Brownian sheet Bq , 
the p-th moment of I(x, y , R ) is equal to 


E[(I(x,y,R)) p ] = P( \\B h (x + P) - B H (y + P) || ^ R~\l ^ j ^ p) dt 1 ■ ■-dt p 

= f P( \B*{x + P) - B»(y + P)\ < iT 1 ,1 < j < p) d dt 1 • • • dt p . 

J[e,l ] Np 

We will bound the previous integral by induction on the parameter p. Hence, let us fix 
fix t 1 ,... , t p_1 e [e, 1]^ and integrate over the variable t p . Note that without any loss of 
generality, we assume that all coordinates of t 1 , ..., t p_1 are distinct. 

The distribution of B^(x + t p ) — Bq {y + t p ) conditionally to B^(x + P) — Bq {y + t p ), 
j e {1,... ,p — 1} is clearly centered and Gaussian. Therefore, 

V{t p ) := P( | B»{x + t p ) - Bq {y + t p )| < R' 1 \ \B*(x + P) - B*(y + P)\, 1 ^ j < p - l) 

< R- 1 • Var (B*(x + t p ) - B*(y + t p ) \ B”(x + P) - B*{y + P), 1< j < p- l) _1/2 

< R- 1 ■ Var«(x + t p ) - B*(y + t p ) \ B”(x + P),B*(y + P),l ^j^p- l)" 1/2 . 


As the reader may expect, we aim to use the anisotropic LND property (2.2) to bound 
the integral V{t p ) dt p . To simplify the former expression, we introduce a collection 

of rectangles (R)i which forms a partition of [e, 1]^ and we split the previous integral 
accordingly. 

More specifically, define for every k e {1,..., TV}, S k = {t{, t 3 k +x k -yk,R k -Xk+yk] 1 < j < 
p— l}. Then, for any index l = (R,... ,Zjv) e {1,... ,3 (jp— 1)}^, let Ii be the TV-dimensional 
rectangle: 


*=n 




+ 



where the elements (s'f : 1 < R < 3(p— 1)) of the set S k are assumed to be increasingly sorted. 
The collection of rectangles (R)i clearly forms a partition of [e, 1]^ (choosing accordingly 
s k and s 3 k p ~ 2 to cover the full square). Consequently, the integration over t p on the domain 
[e, 1]^ can be reduce to a finite sum of integrals on each element R. Thus, let us now set 
l 6 {1,..., 3(p — 1)}^ and observe that for any t p e R, 


Vfce TV}; 


K 


< min | x k + t? 
i^j^p-i 


+ min 

i =si^p- 


, I Vk + t p k 


-k\i 


where z° denotes either x + P or y + P. As a consequence, according to Proposition 3, for 
any t p e R, V(t p ) ^ Cq R~ 1 p(x, y)~ 1 p(s 3 , t p ) _1 , and thus, 

f V(t p ) d dt p ^c 0 R- d p(x,y)- d f p{P ,t p )~ d dt p . 

Jii Jh 


Let us prove the last integral is finite: 

N 


r. r. r TV 'v —a r. TV 

p(s J ,t p )~ d dt p < c j YV| H * f du < c IH|- d ni^| 1/ " / " ldw > 

Jh Jb ( 0 , 1 ) G =1 J Jb(o,i) 


i=i 
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using the simple change of variable vg = u? 1 . Then, switching to spherical coordinates, 
f / 9 (s J , t p ) ~ d dt p < c f r~ d+ ^ t=1 dr f h{ip) dip <+co, 

Jt, Jo Js»-! 


since Yie= 1 V-^ > d and the induced function h is bounded on the sphere S N 1 {1/H^—l > 0 
for every €). Hence, 


f V(tP) d dtP ^c lP N R- d p(x,y)- d , 
J re.ll™ 


de.l] 

and by induction on p, we obtain Inequality (3.1). 


□ 


Remark 3. We may note that the proof of Lemma 3.1 also provides a slighter more general 
inequality. Namely, for any a such that d < a < Tpi 


f 

Jmt 


B*(x + t j ) — Bq (y + t j ) | < R - 1 ,1 < j < p) a dt 1 ■ ■ ■ dt p 

<cg( P !) Jv i2- a M*.y)" ap , 


This extension will be directly used in the proof of Theorem 2. 

The proof of Theorem 1 follows the exact same structure as the ones presented by 
Khoshnevisan et al. [11] and Wu and Xiao [18]. Consequently, we only present the main 
steps, and refer to the former for the technical details which remain the same. 

Proof of Theorem 1. Since B H is Holder continuous with respect to the anisotropic metric 
p : classic results (see for instance [ 20 ] ) on images of fractal sets show that almost surely, 

dim h B h (F + t) < minjd,dim[]F} for all Borel sets E and all t e [0, l] w . 

To obtain the lower bound, we first prove that almost surely, there exists no(w) such that 

Vn ^ noM, Vx, y e [0, 1] N ; I(x, y , 2 n ) sg c 0 n N 2~ nd p(x , y)~ d . 


The previous property is direct application of Borel-Cantelli lemma and the continuity of 
fractional Brownian sheets. We refer to [18] for the details of the arguments. 

Let us now set w e 11, a Borel set F cz [0,1] A , 7 e (0, dim n F) and ye (0, d a 7 ). 
Frostman’s lemma implies the existence of a probability measure p, carrying F and such 
that 

p(S) < ci (diam p S ') 7 for any measurable set S cz [0,1] A '. 

Let v t be the image of p by B H (• + t). Still according to Frostman’s lemma, it is sufficient 
to prove 

CC M d<.Wd») „ 
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to obtain our result. Following the idea of Kaufman [9], we have 
H{dx) p(dy) 


I = 


If || BH( X + t)-BH(y + t)\\v 
V J JJ l[-i,i](R\\B H {x + t) - B H (y+ t)\\)R ri ~ 1 p(dx)p{dy)dR 

V J JJ ![—i,i] (R II B H {x + t)~ B H (y + t)\\)RP~ 1 p(dx)p(dy) d R 


< Id- 


Integrating the previous integral over t e [0,1]^, we thus need to show that 

J := JJ J I(x, y, R)R v ~ 1 dRp(dx)p(dy) 


< oo. 


Let D = {(x,y) e [0, l] 2iv : p(x,y) ^ R 1 } and Ji, J 2 respectively denote the integral J 
over the domains D and D c . Since (p, x p)(D) ^ c\ i? -7 , 

JisSciJ i?" 7+?) - 1 di? < 00 . 

Furthermore, as for any (x,y), I(x,y,R ) < co(w) R~ d p(x,y)~ d , 

Ji < c 0 (o;) \\p(x,y)~ d p(dx)p(dy) f R n ~ d ~ 1 log(R) N dR 
< C 2 (w) JJ p(x, y)~ v log(p(x, y)" 1 )^ p(dx)p(dy) 


< 00 . 


The last two inequalities complete the proof of Theorem 1. 


a 


The second part of this section is devoted to the proof of Theorem 2. As previously, the 
sketch of the latter is highly inspired by the original work of Kaufman [9], and we therefore 
focus on differences compared to the previous results presented by Khoshnevisan et al. [11] 
and Wu and Xiao [18]. 

Proof of Theorem 2. Since dim]] (F) > d, there exists a probability measure p on F such 
that $ r2 n < 00 . To prove that \d(B H (F + t)) > 0, it is sufficient to show that 


.s. I I \V t (u)\ 2 dudt 
J[o,i] w Jr- 1 


< 00 . 


where V t (u) = $ RJV e l ( u - B "( x + t )"> ^(dx) and the exceptional set does not depend on t. 

Let if ^ 0 be a smooth function on such that ip(u) = 1 when 1 < |u| < 2 and if(u) = 0 
outside 1/2 ^ |tt| ^ 5/2. Since S| M | >1 l^( lt )l 2 du is bounded above by 

E 2 n ff ^(2 n B H (x Ft)- 2 n B H (y + t)) p(dx)p(d y), 

„_n o J 


n =0 




Paul Balanga/Image sets of fractional Brownian sheets 


11 


it remains to prove 


V 2 n f (T p(2 n B H (x + t)~ 2 n B H (y + t)) p{dx)p(dy) dt < oo. 

"=o 

For that purpose, we study in the following lemma the component J defined by 
J(x,y,n) := f ip(2 n B H {x + t) — 2 n B H {y + t)) dt. 


n 


Lemma 3.2. There exist two positive constants Cq and f3 such that with probability 1, for 
all n ^ n(uj) and p(x,y) ^ Con N 2~ n , 

J{x,y,n ) (2 + /3)~ n p{x,y)~ d . 


Proof. As observed originally by Kaufman [9], using the reasoning presented in the proof of 
Theorem 1 and a Borel-Cantelli argument, it is sufficient to show the existence of positive 
constants /3, cqj and 03,2 independent of p such that 


E[J(x, y, n) 2p ] ^ 2 + ft)~ 2np p(x, y)~ 2pd . 

Namely, we need to upper bound the following term 

i v 

where t := (i 1 ,..., t 2p ), P e R^. For any n e N, let S n be the following set 


(3.2) 


E 


r 2p 

I ]l\p(2 n B H (x + P )) - P(2 n B H (y + P )) 

Jre.ll 2JV P 7_1 


dt 


2 p AT 


= 


U U{ tS [ £ ’ 1 ] 2iVP : 1^ - 4\ > and 1^ + t k l - 4 - Vt\ > r n Ht 


k= 1 f=l 


where r„ := 03 ^ (71 + 1)2 n . We will begin by studying the former integral on the domain 
S n . It takes the following equivalent form: 


E 


f f Y\exp{i(^ 1 2 n B H {x + P)^2 n B H (y + P))}dfdt 

Js„ Jr 2 ^ j =1 

\ f exp^-^J]Vavfj]^[B H (x + t :, )-2 n B H (y + P)]\ 

JS n J«. 2 r d (. Z J 


2 P 

n *(?)«**■ 


3= 1 


Since 1/2 ^ |£ fe | < 5/2, there exists £q e {1,..., d} such that fff o ^ (2 Vd) 1 . Hence, owing 
the classic sectorial LND property (1.3), 


/ 2 p \ N 

Var( ^ fl[B H (x + P) - 2 n B H (y + P)\ ) ^ c 0 2 2n mini1 1\ - t J e \ 2He , \x t + t\ -P e - ye\ 2Ht \ 

\=i ' i=i 3 * kK 

^ Ci (n + l) 2 . 
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Hence, 

r r ^ 

E FT exp{ , 2 n B H (x + t j ) - 2 n B H {y + t r ))} d£dt < e" 02 ” 2 , 

Js„ JlR 2 *’'* J = i 

where C2 can be chosen sufficiently large up to a modification of 03,0- As a consequence, the 
previous term will clearly be negligible compared to the bound we aim to obtain. 

Let us now consider the second integral over the domain T n := [0, l] 2Arp \S'„, and first 
note that T n can be written as 

2 p n , 

T "=nn { te i £ - i ] 2jvp: 1 < r n Ht } 

jfcetl i=l v 

uIt e [e, l) 2Np : min |o* + t\ - ff 2 - y t \ < . 

I je,2¥=k' )) 

We easily observe that T n is the union of at most (4 n) 2Np sets of the following form: 

Aj = {t e [ £ , l] 2 ^ : \z t + ~ tf k | < r% H ‘, Vfc e {1,..., 2 p}, W e {1,..., N}}, 

where Z£ = 0 or x# — y^. Up to a permutation of indices, the previous set can be written as 
N 

Aj = X{te[e,l] 2p : \z e + t k e - t^ k \ < rl/ H ‘,Vk e {1,... ,2p}}. 

£=1 


As a consequence, using Lemma 3.8 proved in [11], the Lebesgue measure of the previous 
set can be bounded as following: 

A2jv p (Aj) ^2 Np r P /‘ tlir ‘, 


providing a bound on the measure of the full set T n : \ 2 N P (T n ) ^ c^n 2Np r„^ ei Ht ■ Then, 
we divide the integral 


r r ^ i 

I E Yl M2 n B H (x + t j )) -$(2 n B H (y + t j )) 
JT n l j=1 ' 


dt 


into two parts I\ and I 2 , respectively conditioning the former with respect to the events D n 
and D £, where 

D n =\ max \\B H (t + x j )-B H (y + x j )\\>2-^- e)n }. 

U«3«2p M 11 J 

Since i/> is rapidly decreasing, there exists a constant C 5 > 0 which can be chosen as large as 
possible such that 

h^c^n C7P 2 np ^ e=1 H e exp(—C 5 n). 
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Recall that we may assumed that N ^ 2, inducing that > 2, meaning that the 

previous bound is negligible compared to the right end term in Equation (3.2) (/? can be 
chosen small enough). 

Finally, we may conclude the proof by bounding the term I 2 . Let us set a such that 
d < a < Xifli 777 and observe 


f P ( max \\B H (x + r))-B H (y + t j ))\\ 

JT n " / 

= J P( i max \B^{x + t 3 )) - Bq (y + f J ))| ^ 2~^~ e ^ “dt 


i_ — 

A \2N P (T n ) a 


f H max \B*{x + t 3 )) - B*(y + t 3 )) | A 2 -( 1 - £ )")“dt 

Jr £ 1 l2JVp p / 


; [£jl ]2iv P 

using the classic Holder inequality. Based on Remark 3 and the previous estimates, we get 
J P( i max jB H (x + t 3 ))-B H (y + t j ))\\ A 2-( 1 - e)n )dt 
sg cgp c 9"2 " np ( 2d(1_£)+ ( 1_ «) S "= 1 ~k\{x,y)- 2pd 


Then, we may observe that e can be chosen as small as wanted, and particularly, such that 

2 d(l- e) +(l —£) 2 f = 1 £> 2 . 

The combination of the three previous bounds clearly shows the existence of /3 > 0 such 
that 


E[J(a;, y, n) 2p ] sg <^"(2 + P)~ 2np p{x, y )~ 2pd , 

where the constants 03,1 and 03,2 are independent of p e N and neN. Q 

Proof of Theorem 2. To conclude the proof of Theorem 2, we simply observe that the pre¬ 
vious Lemma entails 


00 ft- 00 r- (• 

E 2 " \J(x,y,n)\fj,(dx)fj.{dy) A E 2 ^ 24 ^)”" pO, y)~ d y(dx)y(dy) 

71 = 0 44 


R 2N 


< 00 . 


R 2N 


□ 
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